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1. Introduction 

SMMJ16359+6612 is a sub-millimetre galaxy (SMG) lo- 
cated at z=2.516 that is triply imaged and highly mag- 
nified by the core of a massive foreground cluster, Abcll 
2218 (Kneib et al. 2004a, Knudsen 2004). The three lensed 
images of the background source are highly magnified (by 
factors of 22, 14 & 9) and the maximum image separation 
is ~ 41 arcseconds. All three images are detected by the 
Submillimeter Common-User Bolometric Array (SCUBA) 
at 850 & 450 microns with the brightest image having a 
measured 850 micron flux density of 17 mJy (Kneib et 
al. 2004a). Recently, both Kneib et al. (2004b) and Sheth 
et al. (2004) have observed molecular emission from the 
CO (3-2) line for each of the three images. There are two 
distinct velocity components in this line separated by 280 
kms/s, suggesting a total dynamical mass of the galaxy 
of 1.5 x 10 10 M Q . There is also a spatial offset of 1 arcsec 
between these two components and Kneib et al. (2004b) 
argue, that the source is likely to be a merging system 
with the 2 nuclei separated by 3 kpc. 

At radio wavelengths only the bright tail of the SMG 
population is usually detectable, so observations of highly 
magnified systems such as SMMJ16359+6612, offer a 
unique opportunity to study the radio properties of one 
example of the faint sub-mJy, SMG population. Little is 
known about the nature of these faint SMG, despite the 
fact that they dominate (energetically) the cosmic far- 
infrared background (Knudsen 2004). The intrinsic (un- 
lensed) sub-mm flux of SMMJ16359+6612 is estimated to 
be Ssso — 0.8 mJy (placing it well below the confusion 
limit of normal blank field sub-mm surveys). If one as- 
sumes that the FIR-radio correlation holds for this high-z 
source (e.g. Garrett 2002), the intrinsic flux density of any 
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radio counter-part is estimated to be only a few microJy 
at 8.2 GHz. However, the magnification provided by the 
lens boosts the flux density of even the faintest image in 
SMMJ16359+6612 by a factor of 9, suggesting that ra- 
dio counterparts to the sub-mm sources associated with 
SMMJ16359+6612 should be easily detectable in deep 
Very Large Array (VLA) and Westerbork Synthesis Radio 
Telescope (WSRT) images of the field. We note that singly 
imaged SMG, modestly amplified by foreground clusters, 
have already been detected at radio wavelengths e.g. Smail 
et al. (2000) and Ivison et al. (2001). 

In this paper we present WSRT 1.4 GHz and VLA 8.2 
GHz observations of Abell 2218, specifically the area of 
sky surveyed by the Leiden-SCUBA lens survey (Knudsen 
2004) that includes SMMJ16359+6612. The VLA data are 
in the public domain and were extracted from the on-line 
VLA archive. In section 2 we describe the radio obser- 
vations and data analysis. The results and discussion are 
presented in section 3 and we give our conclusions in sec- 
tion 4. Throughout this paper we assume an Q = 0.3, 
A = 0.7 cosmology with H = 70 km s _1 Mpc -1 . 

2. Radio Observations 

2.1. WSRT Observations 

Observations of Abell 2218 at 1.4 GHz were conducted by 
the WSRT on 5 November 2004. The 12 hour run used 
the full 160 (8 x 20) MHz WSRT observing band and 
employed the default continuum frequency set-up (with 
the eight bands centered between 1.311 and 1.450 GHz). 
Observations were made assuming a coordinate equinox of 
2000. For each of the 8 bands, 64 spectral channels were 
generated (a total of 512 spectral points were obtained 
for the 160 MHz band) and 4 polarization products were 
recorded. Unfortunately, radio telescope "RT7" was being 
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used for single-dish VLBI observations and was unavail- 
able during the entirety of these observations. A short 20- 
minute scan on 3C 48 was used to amplitude calibrate the 
data. A phase reference source (J1642+689) was observed 
for 5 minutes every hour. The data analysis was performed 
using the NRAO AIPS package. The fringe-fitted phase 
solutions from J1642+689 were applied to the A2218 tar- 
get field. For this field, the entire WSRT primary beam 
was imaged using the AIPS task IMAGR. Several bright 
sources were detected in the field of view, and the data 
were self-calibrated using these and other sources in the 
field as a sky model. These outlying sources were later 
subtracted from the data set and self-calibration contin- 
ued using only those sources located in the centre of the 
field. The data were self-calibrated first in phase and later 
in both amplitude and phase using the AIPS task CALIB. 

A uniformly weighted ("robust -1") WSRT image of 
the field convolved with a Gaussian restoring beam of 
13.6 x 12.4 arcseconds (in position angle, PA = -0.21°) 
is shown in Figure 1. At 1.4 GHz there is a great deal 
of extended emission in the vicinity of the A2218 cluster 
core. In order to enhance the contrast of the WSRT im- 
age, spacings shorter than 800A were not used to form the 
image presented in Figure 1. The image reaches a 1 a rms 
noise level of ~ 15 /xJy beam -1 . 



2.2. VLA observations 

We searched the VLA archive for observations associated 
with the central area of Abell 2218. Several data sets 
were identified that included short, snapshot observations 
of this field. Only one full-track 8.2 GHz observation of 
A2218 was identified, AR045. These X-band observations 
were made on 1999 June 18 UT with the VLA in its D (3 
km) configuration. Data were acquired in dual circular po- 
larizations with bandwidths of 25 MHz and at center fre- 
quencies of 8.1732, and 8.2732 GHz. Each of these 25 MHz 
IFs was further divided into seven spectral channels with a 
width of 3.125 MHz. Observations were made assuming a 
coordinate equinox of 2000. The absolute flux density scale 
was set by observations of 3C 286. The data calibration 
was made in the standard way with the phase variations 
during the observations calibrated via short observations 
of the VLA calibrator source 1642+689. Unfortunately 
observations of 1642+689 were only made 3 times over 
the course of the 24 hour observations. The calibrated 
A2218 data were clipped in order to remove a few high 
points before beginning the self-calibration process. The 
self-calibration process followed the same outline as that 
described for the WSRT observations. 

A naturally weighted VLA image of the field convolved 
with a Gaussian restoring beam of 11.2 x 10.4 arcseconds 
(in position angle, PA = —72°) is shown in Figure 2. The 
image has an r.m.s. noise level of 6 microJy beam. 
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Fig. 1. The WSRT 1.4 GHz contour map superimposed 
upon an HST F702W image of the region of sky associated 
with SMM J16359+6612. Radio emission from all three 
lensed images (SMM-A, SMM-B and SMM-C) is detected 
by the WSRT. The boxes represent the SCUBA source po- 
sitions. Radio emission from three other SCUBA sources 
in the field (including the highly magnified arclet #289, 
SMM J16359+66118) are also boxed. Contours are drawn 
at -3,3,5,10,20 & 40 times the l-er noise level of 15 mi- 
cro Jy/beam. 

3. Results & Discussion 

Figures 1 and 2 present the uniformly weighted WSRT 
and naturally weighted VLA contour maps of part of the 
A2218 field, superimposed upon an HST F702W image 
(Kncib ct al. 1996) of the region of sky that includes 
the lensed sub-mm galaxy SMM J16359+6612 (Kneib et 
al. 2004a). The positions of the nine sources detected by 
SCUBA at 850 micron (Knudsen 2004) arc also identi- 
fied as square boxes in these figures. Following Kneib 
et al. (2004a) we identify the lensed images as SMM- 
A (SMMJ16359+6612.6), SMM-B (SMMJ16359+6612.4) 
and SMM-C (SMMJ16358+6612.1). Both the WSRT and 
VLA observations detect radio emission associated with 
all three lensed images. The faintest SCUBA detection 
(SMM-C) is detected at the 4-sigma level by the VLA at 
8.2 GHz and at the 7-sigma level by the WSRT at 1.4 
GHz. 

The AIPS task IMFIT was used to fit Gaussian com- 
ponents to the WSRT and VLA radio sources. Details of 
the radio source's total flux density (St), peak flux den- 
sity (Spfc), position and (where appropriate) deconvolved 
Gaussian sizes (major axis, minor axis and position an- 
gle) are detailed in Table 1, together with their associated 
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Table 1. Details of the WSRT and VLA radio sources associated with SMM J16359+6612 (SMM-A, B, C) and 
SMM J16359+66118 (arclet #289). 
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WSRT J163550+661205 (SMM-C) 


51.1 ±0.1 


12, 08.0 ± 1.0 


110 ± 17 


110 ± 17 








WSRT J163553+661226 (SMM-B) 


53.9 ±0.1 


12, 26.4 ±0.6 


391 ± 43 


217 ± 16 


15 ±2 


9±2 


165 ± 12 


WSRT J163555+661237 (SMM-A) 


55.6 ±0.2 


12, 37.7 ± 1.0 


99 ± 17 


99 ± 17 








VLA J163550+661206 (SMM-C) 


50.5 ±0.2 


12, 06.8 ±0.1 


24 ± 10 


24 ±6 








VLA J163554+661223 (SMM-B) 


54.5 ±0.1 


12, 25.3 ± 1.1 


124 ±18 


51 ±6 


16 ±3 


10 ±3 


39 ± 12 


VLA J163555+661236 (SMM-A) 


55.7 ±0.1 


12, 37.1 ±0.5 


28 ±6 


28 ±6 
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Fig. 2. The VLA 8.2 GHz contour map superimposed 
upon an HST F702W image of the region of sky as- 
sociated with SMM J16359±6612F702W. Radio emis- 
sion from all three lensed images is detected by the 
WSRT. The boxes represent the SCUBA source posi- 
tions. Radio emission from the highly magnified arclet 
(#289), SMM J16359±66118 is also detected. Contours 
are drawn at -3,3,4,7 & 10 times the 1-a noise level of 6 
micro Jy/beam. 



formal errors. It was not possible to robustly constrain 
the deconvolved size of the Gaussian component associ- 
ated with the fainter images SMM-A and SMM-C in the 
VLA and WSRT images. For SMM-A and SMM-C the 
Guassian component size was thus fixed to the respective 
restoring beam sizes. For both the VLA and WSRT im- 
ages, two Gaussians were fitted to the SMM-A, B region 
simultaneously. 



The SCUBA, VLA and WSRT observations have com- 
parable resolution. The error in the SCUBA 850 micron 
positions (including confusion effects) is expected to be 
~ 4 arcseconds (Knudsen 2004). The radio positions are 
expected to be better than this, with an accuracy of 2 arc- 
second or better. Slightly larger errors are possible for the 
radio positions derived by the WSRT, due to mild confu- 
sion in the field. The radio and sub-mm positions of SMM- 
A,B,C are all consistent within the errors - the largest 
deviation is an offset of 3 arcseconds (the VL A/SCUBA 
comparison in the case of the faintest image SMM-A). 
There are no obvious systematic offsets in the sub-mm 
and radio positions. 

At 8.2 GHz the total flux density of SMM-B is 124 
microJy, somewhat larger than our original expectations 
(see section 1). By co- adding the total flux density of all 
three images at 1.4 and 8.2 GHz, we derive a spectral in- 
dex of the source a = —0.7, similar to the steep values 
measured by Richards (2000) for star forming galaxies in 
the Hubble Deep Field North. Kneib ct al. (2004b) have 
shown that the overall SED of SMMJ16359±6612 is sim- 
ilar to Arp 220. By adopting an Arp 220 SED, we esti- 
mate the source's k-corrected radio luminosity at 1.4 GHz 
to be ~ 4 x 10 23 Watts/Hz. Using the relation between 
Star formation rate (SFR) and radio luminosity (Condon 
1992), we derive an intrinsic star formation rate (SFR) for 
SMM J16359±6612 of 500 M Q yr" 1 . This is in good agree- 
ment with the value obtained by Kneib et al. (2004) from 
the FIR luminosity. As noted by Kneib et al. the uncor- 
rected Ha SFR estimate is only 11 M yr -1 , suggesting 
that this galaxy is highly obscured by dust. 

The VLA 8.2 GHz observations best resolve the lensed 
images. The flux density ratio of the images at 8.2 GHz 
(and their formal error) are ~ 0.23±0.06 and ~ 0.19±0.08 
for SMM-A/SMM-B and SMM-C/SMM-B respectively. 
Similarly, the WSRT 1.4 GHz flux image density ratios 
are 0.25 ± 0.05 and 0.28 ± 0.06 respectively. The radio 
flux density measurements at 1.4 and 8.2 GHz are thus 
consistent with the source being gravitationally lensed. In 
comparison, the SCUBA observations give higher values of 
0.64 and 0.53, the IRAM CO(3-2) line intensities (Kneib 
et al. 2004b) give 0.67 and 0.63, and the K-band observa- 
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tions (Kneib et al. 2004a) give 0.53 and 0.4. Lens models 
of this system (based on Kneib et al. 1996), predict flux 
density ratios of 0.63 and 0.4. Discrepancies in the flux 
density ratio at different wavelengths might be explained 
if the various emission regions are not co-located or have 
very different size-scales over which the magnification may 
change. We also note that the flux density ratio of SMM- 
A, B and C between 850 micron and 1.4 GHz are 111 ±22, 
43±7 and 82±16. By comparison the 850 micron/1.4 GHz 
flux density ratio of another SMG also located at z ~ 2.5, 
SMM J1401 1+0252 (Ivison et al. 2000) is 127 ±37. It may 
be that in the case of SMM J16359+66I2, the discordant 
flux ratios at both 1.4 and 8.2 GHz can be explained by 
enhanced radio emission associated with SMM-B, perhaps 
from another radio source also located in this region of the 
field. Higher resolution radio observations are required in 
order to settle this question. 

As shown in Table 1, the brightest image SMM-B is 
resolved by the VLA 8.2 GHz observations. The measured 
PA of the major axis is ~ 39 degrees, in good agreement 
with the position angle of the major axis of the CO (3-2) 
measurements of the same component, ~ PA = 30 de- 
grees. Both measurements are consistent with the overall 
extension of the associated arc-like HST images (Kneib et 
al. 2004a). The FWHM of the major axis of the SMM- 
B (see Table 1) does not exceed <~ 17 arcseconds. With 
a magnification of ~ 22 for this image, the intrinsic size 
must be on the arcsecond or sub-arcsecond scale - consis- 
tent with the sub-galactic sizes measured for radio sources 
in the HDF-N. 

The WSRT observations also detect other radio coun- 
terparts to other sub-mm sources detected by SCUBA 
in this field (Knudsen 2004). In particular, we detect ra- 
dio counter-parts to SMMJ16357±66117 (4cr detection), 
SMM J16361+66126 & SMM J16359+66118. The latter 
source (see Table 1) is previously identified by Kneib et al. 
(1996) as a z=1.034 singly-imaged lensed arc (also known 
as arc #289, see Swinbank et al. 2003 and references 
therein) with a magnification of ~ 7. The source is also de- 
tected by the 8.2 GHz VLA observations presented here, 
and by the ISOCAM at 15 micron (Barvainis, Antonucci 
& Helou 1999). In both the WSRT and VLA images, the 
measured size of the arclet is less than the restoring beam, 
in particular the minor axes of the fitted Gaussians arc un- 
resolved. The position angle of the major axes are consis- 
tent with the optical extension of this system. This source 
also has a steep spectral index, a ~ 0.9. The VLA 8.2 GHz 
observations do not detect the other SCUBA sources - 
SMMJ16357±66117 & SMM J16361±66126. These are 
fainter than SMM J16359+66118, and must have fairly 
steep spectral indices (a > 0.5). Further details of the ra- 
dio counterparts to these SCUBA sources will be reported 
elsewhere. 

4. Conclusions 

We have detected with the WSRT at 1.4 GHz and VLA at 
8.2 GHz, radio emission associated with the triply lensed 



SMG, SMMJ16359±6612. This is the first time that ra- 
dio emission has been detected in a multiply imaged SMG 
system lensed by a foreground cluster. The maximum im- 
age separation is ~ 41 arcseconds, much larger than any 
other lens system detected in the radio. This is also the 
first time that multiply lensed radio emission has been de- 
tected from a star forming galaxy — all previous multiply 
imaged radio lensed systems are associated with radio- 
loud AGN. The properties of the three radio sources are 
largely consistent with the gravitational lensing hypoth- 
esis. In addition to detecting SMMJ16359+6612, we also 
detect radio emission coincident with SMMJ16359+66118 
a singly imaged arclet (#289). The sources are only de- 
tectable by current radio instruments due to the high mag- 
nification factor provided by the lens (^ 45 in the case of 
SMMJ16359±6612). The intrinsic total flux density of the 
radio source at 1.4 and 8.2 GHz is ~ 14 and ~ 3 microJy. 

Follow-up radio observations of SMMJ16359±6612 
with much better angular resolution than those presented 
here are warranted. Since sub-arcsecond imaging is rou- 
tinely possible at radio wavelengths, it should be possi- 
ble to fully capitalise on A2218 as a "natural telescope". 
A comparison with the detailed HST images of this sys- 
tem will be particularly interesting. Measurements of the 
size and a study of the morphology of the radio emit- 
ting region should also constrain the extent of the star 
formation in this galaxy. Upgraded radio telescopes (such 
as the EVLA and e-MERLIN) should be able to detect 
many more highly magnified lens systems lying behind 
massive foreground clusters. In very deep integrations, the 
c-MERLIN and the EVLA telescopes now under develop- 
ment, may be able to detect radio sources associated with 
individual SNe (or SNR), assuming these have luminosi- 
ties similar to those detected in local mergers, such as Arp 
220. Radio observations of SMMJ16359±6612 and other 
highly magnified SMG using existing and upgraded ra- 
dio telescope facilities, will be important in determining 
the general radio properties (e.g. angular size) of the faint 
SMG population, providing essential input to the techni- 
cal specification of next generation instruments, such as 
the Square Kilometre Array. 
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